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Within many-body perturbation theory we apply vertex corrections to various closed-shell atoms and to jel- 
lium, using a local approximation for the vertex consistent with starting the many-body perturbation theory from 
a DFT-LDA Green's function. The vertex appears in two places - in the screened Coulomb interaction, W, and 
in the self-energy, E - and we obtain a systematic discrimination of these two effects by turning the vertex in E 
on and off. We also make comparisons to standard GW results within the usual random-phase approximation 
(RPA), which omits the vertex from both. When a vertex is included for closed-shell atoms, both ground-state 
and excited-state properties demonstrate little improvement over standard GW. For jellium we observe marked 
improvement in the quasiparticle band width when the vertex is included only in W, whereas turning on the 
vertex in E leads to an unphysical quasiparticle dispersion and work function. A simple analysis suggests why 
implementation of the vertex only in W is a valid way to improve quasiparticle energy calculations, while the 
vertex in E is unphysical, and points the way to development of improved vertices for ab initio electronic struc- 
ture calculations. 

PACS numbers: 71.45.Gm, 31.25.Eb, 31.25Jf, 71.10.Ca 



I. INTRODUCTION 

Many-body perturbation theory (MBPT) is a leading 
method for computing excited-state electronic properties in 
solid-state physic s 1 ' 2 ' 3 . Within many-body perturbation the- 
ory, Hedin's GW method 4 is the most widely used approxima- 
tion for the self-energy, E. The exact one-body Green's func- 
tion, G (which contains information about ground and excited- 
state properties of the system) can be written, using a Dyson 
equation, in terms of a suitable Green's function of a "zeroth- 
order" system of non-interacting electrons, Go (constructed 
from that system's one-particle wavefunctions and energies), 
and the self-energy operator E. The approximation is defined 
by the choice of zeroth-order system and by the expression 
(typically a diagrammatic expansion in terms of Go) used to 
approximate E. The self energy, E, contains all the informa- 
tion of many-body interactions in the system and can be ob- 
tained by using Hedin's set of coupled equations: 

E(12) = i [ W(l + 3)G(14)r(42;3)<i(34) (1) 



W(12)=v(12)+ /W(13)P(34)v(42)J(34) (2) 



P(12) = -i / G(23)G(42)r(34; l)c/(34) (3) 
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r(12;3) =5(12)5(13) 

+ /|||^yG(46)G(75)r(67;3) £ /(4567), (4) 

and the Dyson equation, where P is the polarizability, W the 
screened and v the unscreened Coulomb interaction and T the 
vertex function. The notation 1 = (xi,<7i,fi) is used to denote 
space, spin and time variables and the integral sign stands 
for summation or integration of all of these where appropri- 
ate. (1 + denotes t \ + TJ where 77 is a positive infinitesimal 
in the time argument). Atomic units are used in all equa- 
tions throughout this paper. These are four coupled integro- 
differential equations where the most complicated term is the 
vertex, F, which contains a functional derivative and hence, 
in general, cannot be evaluated numerically. The vertex is the 
usual target of simplification for an approximate scheme. 

The widely used GW approximation is derived with the 
Hartree method as a starting point, and hence has a rig- 
orous foundation only when started from a non-interacting 
Green's function, Go, made from eigenstates of the Hartree 
Hamiltonian. This is because the initial self-energy, Eo = 
and the vertex function is correspondingly set to T(12;3) = 
5(12)5(13) since 5E(12)/5G(45) = 0. 

Solving Hedin's equations with the vertex fixed in this ex- 
pression yields the so-called self-consistent GW approxima- 
tion. In this approach, the self-energy operator is formed from 
a product of a Green's function and a screened Coulomb inter- 
action, where the Green's function used is consistent with that 
returned by Dyson's equation. Since the self-energy depends 
on G, this procedure should be carried out self-consistently, 
beginning with G = Gq. 

In practice, it is customary to use the first iteration only, 
often called GoW RPA , to approximate the self-energy opera- 
tor. Here, W RPA is perhaps the simplest possible screened 
interaction, which involves an infinite geometric series over 
non-interacting electron-hole pair excitations as in the usual 
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definition of the RPA. It is important to make this one itera- 
tion as accurate as possible, so an initial Go calculated using 
Kohn-Sham density-functional theory in the local-density ap- 
proximation (DFT-LDA) is normally used. This choice of Go 
generally produces much more accurate results for quasipar- 
ticle energies (the correct electron addition and removal en- 
ergies, in contrast to the DFT-LDA eigenvalues 5 ). However, 
because this choice of Go corresponds to a non-zero Eq, there 
is no longer a theoretical justification for the usual practise of 
setting the vertex to a product of delta functions and different 
choices for the exchange-correlation functional may lead to 
different Green's functions^. 

Using the static exchange-correlation kernel, K xc , (which is 
the functional derivative of the DFT exchange-correlation po- 
tential, y xc , with respect to density, n) Del Sole et a/. 8 demon- 
strated how GoW RPA may be modified with a vertex function 
to make E consistent with the DFT-LDA starting point. They 
added the contribution of the vertex into both the self-energy, 
E, and the polarization, P. The result is a self-energy of the 
form GoWbr LDA9 . The GoW LDA approximation is obtained 
when the vertex function is included in P only. As com- 
mented by Hybertsen and Louie 10 and Del Sole et al. , both 
these results take the form of GW, but with W representing the 
Coulomb interaction screened by respectively the test-charge- 
electron dielectric function and the test-charge-test-charge di- 
electric function, in each with electronic exchange and cor- 
relation included through the time-dependent adiabatic LDA 
(TDLDA). 

Del Sole et al. found that G Wor LDA yields final results 
almost equal to those of GoW RPA for the band gap of crys- 
talline silicon and that the equivalent results from GoW LDA 
were worse when compared to GoW RPA . It should perhaps be 
mentioned that the inclusion of other types of vertex correc- 
tions have been studied before as well, most notably correc- 
tions based on various approximations of a second iteration 
of Hedins equations, starting with GoW RPA1 ^. However, 
these have usually been applied with initial Kohn-Sham (KS) 
Green's functions, which are still not theoretically consistent 
with that starting point. The correct theoretical treatment of 
a second-iteration vertex from KS Green's functions is quite 
complicated and still absent in the literature. 

The purpose of the present work is to make a systematic 
study, for both localized and extended systems, of a simple 
ab initio vertex correction whose form is determined by the 
starting approximation for the self-energy (Eo = V xc for DFT- 
LDA). Related vertex corrections, including others derived 
from K xc , have been investigated in earlier work. For exam- 
ple, Northrup et al. 13 used LDA bulk calculations as a starting 
point and a plasmon-pole calculation of the response func- 
tion in conjunction with a GoW LDA -like vertex correction in 
the screened interaction. They found a narrowing of the band 
widths of Na, Li and Ki 4 in agreement with the experiments 
of Jensen and Plummer 15 who had noted that the experimental 
band width was significantly narrowed (~ 23 %) compared to 
the free-electron result. Hedin's GoW RPA4 calculations only 
gave a narrowing of about 10% for an homogeneous electron 
gas of the same mean density, indicating a large impact of 
further many-body effects. This led to additional experimen- 



tal and theoretical investigation s 1 6 i 1 7 i 1 8 i 19 ' 20 - 2 1 but the issue re- 
mains controversia l 22 i 23 ' 24 i 25 . 

For individual atoms, GW quasiparticle properties have 
been investigated previously by Shirley and Martin 26 (includ- 
ing an exchange-only vertex) and, more recently, total energy 
studies on atoms and molecules using the variational function- 
al of Luttinger and Ward 27 have been performed by Dahlen 
et aL 2 ^ 29 , Stan et al^- and Verdonck et al.—. These stud- 
ies have shown that GoWq^ in general gives quasiparticle 
properties which are much improved over DFT and Hartree- 
Fock methods and that, when calculated self-consistently, GW 
also provides reasonably good total energies for atoms (with 
differences versus highly accurate reference methods being 
on the order of tens of mHa per electron). To its merit, 
self-consistent GW is also a conserving approximation in the 
Baym-Kadanoff 32 sense. However, non-self consistent total 
energies in GoW RPA are noticeably less accurate. Conversely, 
the good agreement between the quasiparticle energies and 
experiment is destroyed when performing self-consistent cal- 
culations. 

The answer to why this happens must, by definition, lie with 
the only approximated quantity, the vertex correction. This 
study is meant to address the need for a precise (including a 
full treatment in frequency) comparative study of the vertex 
corrections proposed by Del Sole et al. for localized and ex- 
tended systems within GoWq. 



II. METHOD 

Hybertsen and Louisi^ comment that it is possible to 
start a GW calculation from an initial self-energy, Eo(12) = 
5(12)V XC (1). This approach gives a theoretical basis for be- 
ginning a GoWo calculation from DFT-LDA orbitals. Adopt- 
ing this idea, we see from Eq. (0]l that the second term is now 
non-zero, unlike in the GW approximation. Since the electron 
density is n(\) = — iG(l 1 + ) then, 



5E(12) _ <5E(12) 8n(4) 
<5G(45) 



<Sn(4) 5G(45) 
..5V X L C DA (1) 



5«(4) 



(5) 
(6) 



where delta-functions are to be understood in all other vari- 
ables. In an appendix Del Sole et al} show how to add this ap- 
proximate vertex to both W and E, and into W only, by form- 
ing two different effective Ws. Our method follows that of Del 
Sole et al.— by modifying the dielectric function, e from its 
form in the RPA. The screened Coulomb interaction in MBPT 
is written as 



W 



(7) 



where g -1 is the inverse dielectric function. We use the full 
polarization without recourse to plasmon-pole models. The 
random phase approximation (RPA) dielectric function is 



(8) 
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Del Sole et al. show that adding the form of the vertex from 
Eq. (|6]l into both E and W modifies the RPA dielectric function 
to, 



~e = l-(v + K xc ) X c 



(9) 



which leads to the introduction of an effective screened 
Coulomb interaction W . This is trivial to implement into a 
GW computer code as it requires a simple matrix addition, 
once K xc is calculated. The result of this modification is 
that W contains not only the screened Coulomb interaction 
but also an exchange-correlation potential. We shall refer to 
this method as GoWqT lda as we have added the correct DFT- 
LDA vertex to the GW method, hence the method is a one- 
iteration GWT (GqWoTq) calculation beginning with a DFT- 
LDA Green's function. 

An alternative choice for the effective dielectric function, 
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corresponds to adding K xc into W only. We term this method 
GqWq DA as the LDA vertex is inserted into the screened 
Coulomb interaction, W, only. This is equivalent to the one- 
iteration GW approximation, Go Wo but with W calculated us- 
ing the adiabatic LDA. 

The vertices presented are thus dynamical, i.e. frequency 
dependent, due to the inclusion of % , and must include the 
excitonic effects of the corresponding jellium due to the ap- 
pearance of K xc . Another way of looking at it is that this cor- 
responds to a treatment beyond GqWq where at the level of 
the vertex corrections, the system is modelled by the homoge- 
nous electron gas. It is not likely, however, that these methods 
would be able to capture any satellite structure beyond that 
provided by K xc as the calculations are non-self-consistent. 



III. COMPUTATIONAL APPROACH 

The quasiparticle energies, £,■ and wavefunctions, l//,-, are 
formally the solution of the quasiparticle equation, 

{-IV 2 + V ext (r)+V H (r)} y/ ( -(r)+ /E(r,r';£ ; )y/,(rVr' 

(11) 



where Vext and Vh are the external and Hartree potential, re- 
spectively. 

In the case of a spherically symmetric system it is sufficient 
to describe all non-local operators in the GW formalism by 
two radial coordinates and one angular coordinate, 0, that de- 
notes the angle between the vectors r and r'. The self-energy, 
E, then assumes the much simpler form, 

E(r/,0;e)=£[E,(r/;e)]P z (cose), (12) 

1=0 

where P/(cos 9) is a Legendre polynomial of order /. 

The Legendre expansion coefficients of the self-energy, E, 
are calculated directly, thereby circumventing the need for a 



numerical treatment of the angular dependence. We use a 
real-space and imaginary time representation^ to calculate 
the self-energy from the non-interacting Green's function Go- 
The self-energy on the real frequency axis, required for solv- 
ing the quasiparticle equation, is obtained by means of an- 
alytic continuation^. The current implementation has been 
successfully applied to jellium clusters 34 and light atoms 7 ' 35 . 

To obtain the quasiparticle energies and wavefunctions the 
quasiparticle equation ( fTTT i is fully diagonalized in the basis of 
the single particle orbitals of the non-interacting Kohn-Sham 
system. For localized systems the quasiparticle wavefunc- 
tions can differ noticeably from the wavefunctions of the non- 
interacting system or in certain cases even have a completely 
different character, as was demonstrated for image states in 
small metal clusters^ 4 -. 

Ground-state total energies were calculated using the 
Galitskii-Migdal formula 36 transformed to an integral equa- 
tion over imaginary frequency. This avoids analytic continu- 
ation of the self-energy, which would be unreliable for large 
frequencies. 

For jellium, the homogenous electron gas, we solve Hedin's 
equations in wavevector and real-frequency space. This 
avoids analytic continuation and enables accurate and easy ex- 
traction of spectral properties. Again, we do not use plasmon- 
pole models, but the full frequency-dependent polarization. 



IV. TOTAL ENERGIES 



Method 


He 


Be 


Ne 


HF 
DFT-LDA 

GoW$* A 

G iv LDA 
G w r LDA 

VMC 
DMC 
CI 


-1.4304" 

-1.4171 

-1.4117(5) 

-1.4120(2) 

-1.3912(2) 

-1.45176° 

-1.45186" 

-1.45189'' 


-3.6433 s 

-3.6110 

-3.5905(9) 

-3.590(1) 

-3.573(1) 

-3.66670'' 

-3.66682 c 

-3.66684 rf 


-12.8547" 

-12.8183 

-12.777(1) 

-12.775(15) 

-12.745(10) 

-12.891(5)" 

-12.89231" 

-12.89370"' 



" See reference 
b See reference 
c See reference 
d See reference 



TABLE I: Total energy data (Hartrees/electron). See Fig.Q] 
— A comparison of various methods for total energy calcula- 
tions. Hartree-Fock, Density-Functional Theory, one-iteration GW 
(GoWq RPA ), the two approximate vertex GWs, variational Monte 
Carlo, diffusion Monte Carlo and configuration interaction. (CI usu- 
ally yields the most accurate estimate of the ground-state energies for 
localized systems.) 

The MBPT total energy results are compared against con- 
figuration interaction (CI) and quantum Monte Carlo methods 
(variational Monte Carlo (VMC) and diffusion Monte Carlo 
(DMC)). The CI and QMC family of methods usually yield 
the most accurate estimates of ground-state energies and are 
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FIG. 1 : Total energies of atoms. - We compare a number of GW- 
based approaches to Hartree-Fock 37 ; DFT with an LDA exchange- 
correlation fucntional consistent with the A^c used; quantum Monte 
Carlo 3 ! (VMC and DMC) and CM. The dotted line is the CI value 
and is there to guide the eye. In all cases GoW r LDA behaved poorly 
in comparison to GoW ( RPA , whereas GqWq' makes no improve- 
ment to GqWq . The MBPT methods are not as accurate as the 
computationally cheaper mean-field calculations but GqWq and 
GqW ( \- da are the better of the three MBPT methods. 





2 


3 


4 


5 


G < PA 
G W LDA 
G W4,r LDA 
GW a 
QMC (DMC) b 


-0.2826(3) 
-0.2857(4) 
-0.2525(2) 
-0.2727(5) 
-0.2742(1) 


-0.1967(1) 
-0.2002(2) 
-0.1678(1) 

n/a 
-0.1902(1) 


-0.1522(1) 
-0.1560(1) 
-0.1241(1) 
-0.1450(5) 
-0.1464(1) 


-0.1247(1) 
-0.1288(1) 
-0.0972(1) 
-0.1185(5) 
-0.1202(1) 



" See reference 
b See reference 
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TABLE II: e xc (Ha) — The exchange correlation energy for jel- 

|£f + £ X c, where 



lium. The total energy per particle is given by e : 

£ F = jfcp = j (w^) anl ^ a ~ iwa ) 3 • ^ e ener gi es under heading 
GW are from self consistent calculations by Garcfa-Gonzalez and 
Godby 41 for reference. GqWq^ is lower than the QMC energy by 



- 5% on average. G W LDA is ' 



- 6 % lower and G Wbr LDA is ' 



10% 



too high. The DMC values are evaluated by Perdew and Zunger 
parametrization of Ceperly and Alder's DMC calculations. 



A2 



three cases. It is known that there is a large self-interaction 
error in the LDA, especially noticeable in smaller atoms. 
Hartree-Fock, which is self-interaction free by construction, 
is more accurate than DFT-LDA. Hence the self-interaction 
error is introduced via the LDA orbitals into the Green's func- 
tion, Gq DA , which gives rise to the GqW^^ total energy's 
consistent poor agreement with CI. (By way of illustration, 
using a Go from the superior KLI 43 , an optimized effective 
potential that is formally free of self-interaction error for a 
two-electron system, greatly improves the DFT and GW to- 
tal energies. The Gq^Wq^ results for He, Be and Ne are 
-1.4550(3), -3.6780(2) and -12.843(1) respectively.) 

For all three atoms the vertex in W alone (GoW LDA ) shows 
little difference to G W RPA (TableU and Fig.Q}, whereas 
GoWoT LDA raises the total energy with respect to GqWq . 
This change is due not to the LDA self-interaction but to 
the nature of the vertex. The result of adding the LDA ver- 
tex to G^Wq mirrors that of adding it to Gq DA W . There 
is an increase of the total energy when the vertex is applied 
in both W and £ (G^Wbr LDA ) but the vertex in W only, 
(G^W^), results in a similar total energy to G^W^. 
(The G^WqT^ and G^W^ DA for He are -1.4235(10) 
and —1.4475(5) respectively.) 

In jellium the trend is the same for all densities in the re- 
gion from r s = 2 to 5 (r s is the density parameter, where 

1 /3 

rs= ( ? l_) i / J and n is the electron density in atomic units) as 
can be seen in TableHII GoW LDA lowers the total energy of jel- 
lium slightly as compared to GoWq~ pa and GoWqT lda makes 
the energy too high. GqW^^ is on average ^5% lower than 



LDA 



the QMC result. G W r LUA is - 10% too high and G W L 
~ 6% lower than the QMC result. 

For jellium, neither method leads to a result more accurate 
than GoWq RPA . It is apparent, however, that the vertex added 
solely in the polarization has the minor effect of lowering the 
total energy. When the vertex is subsequently added into the 
self-energy there is a major positive shift in the total energy 
as seen in the atomic results as well. Self-consistent GW 
calculations 41,44 for jellium show that the self-consistent total 
energy is about 4 — 5% higher than the GqW^ a ones in the 
range of r s = 2 to 5 and the essentially exact QMC energies 
are about 0.5 — 1% lower than the self-consistent GW values. 
Assuming to a first approximation that the vertex corrections 
are independent and additive corrections to self-consistency, 
this would indicate that the GoWqT lda energies would still be 
much too high, but the GoW LDA energies might end up very 
close to the QMC results if self consistency is achieved, since 
they lower the GoWq PA energies to roughly the same extent 
as the difference between QMC and self-consistent GW ener- 
gies. 



V. QUASIPARTICLE ENERGIES 



variationally bound, meaning that the lowest energy is the 
most accurate. 

The GqWq result with Go constructed from DFT-LDA 
eigenstates, (GoWq RPA ) is in poor agreement with CI in all 



The quasiparticle energy corresponding to the first ioniza- 
tion energy 63 is presented for helium, beryllium and neon 
in Fig. |2] The MBPT methods are consistently more ac- 
curate than DFT-LDA Kohn-Sham eigenvalues. However, 
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Method 


He 


Be 


Ne 


DFT-LDA 

G ( < DA 

G w r LDA 

CI 

Experiment 


-15.4877 

—74 70(41 

-24.05(5) 

-22.5(1) 

-24.5930° 

-24.587* 


-5.5909 

—9.24(2) 

-9.25(3) 

-7.56(6) 

-9.3226" 

-9.3227 c 


-13.503 
—70 55(101 

iU.JJ 1 1UI 

-19.48(10) 
-18.85(5) 
-21.6034" 
-21.5645°' 




a See reference 
b See references J45B46 
c See reference |47 
d See reference 



TABLE III: First ionisation energy (eV) - A comparison of various 
methods for quasiparticle energy calculations: DFT-LDA, G W RPA 
and the two approximate GWs. CI denotes the ionisation potential 
calculated from the difference in CI total energies and Experiment 
the measured value. 



again GoWq DA is roughly equivalent to GoW RPA for helium 
and beryllium and in all cases GoWbr LDA causes an increase 
in quasiparticle energy, in agreement with Del Sole et a/A 
In general, GqWqT lda worsens QP energies as compared to 
G < PA 

For jellium, different quantities are accessible at differ- 
ent stages of the iteration of Hedin's equations. The pair- 
correlation function g(r) for example, can be obtained from 
the (isotropic) inverse dielectric function, e^ 1 (k,co), by the 
integration 



2rk\ Jo 



dkksm(kr)[S(k)-l], (13) 



where the static structure factor, 



S(k) = -^r d ^3[e-\k,a>)]. 
5 Jo n 



(14) 



g(r) is shown in Fig. [3] for r s — 1.0. The RPA displays the 
well known failure to be positive definite for r s > 0.78. This 
is remedied by the local vertex, but the result appears to be an 
overcorrection (note that GoW LDA and GoWqT lda are equiv- 
alent at this stage since E has not yet been calculated). 

The tendency of GoWoT LDA to overshoot - the reason for 
which, we will defer to the closing discussions - is apparent 
in all subsequent results. Once E has been calculated, the QP 
dispersion can be extracted. Presented in Fig.|4]is the real part 
of the self-energy evaluated at the self-consistent eigenvalues, 
i.e. the correction Re [E k (£k)] to the quasiparticle dispersion 
as found by the formula 



£ k = e£+/?e[E k (£ k )] : 



(15) 



where is the non-interacting dispersion. Care has been 
taken to align the Fermi energy of non-interacting and in- 
teracting systems so that the Dyson equation is consistent 36 
and all quantities are calculated in real frequency. The self- 
consistent quasiparticle energy should be used when one has 
a self-consistent E, but for a GqWq calculation there is still 
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FIG. 2: QP energies of atoms - The first ionization energy is plot- 
ted. We compare the GW based approaches to DFT-LDA; the 
experimenta l 45 ! 46 ! 47 answer (Expt.) and the values calculated from 
the differences in Ctf calculations (CI). The helium, beryllium and 
neon values are plotted with circles, diamonds and squares respec- 
tively. The dotted lines go through the CI value and are there to 
guide the eye. All of the GW calculations are more accurate than 
the mean-field method. For helium and beryllium the GqW rpa and 
GoWq jDA methods give similar results. In all cases GqWoT is 
shifted to a higher energy that G W RPA . (This shift for GoW r LDA 
was also found by Del Sole et a/. 8 and Fleszar and Hanke 49 .) 



controversy about whether the self-consistent eigenvalues £ k 
or the zeroth-order eigenvalues e k ' are best used as the argu- 
ment of E k in equation (TT3t 1,5Q - The self-consistent approach 
was chosen in this paper. 

The difference between the quasiparticle energies at k = kp 
and k = is known as the band width, which therefore takes 
the form of the free-electron value (fep/2) corrected by the 
change in Fig. |4]between k~kp and k = 0. This band width is 
shown in Table|IV|and Fig. [5] It consistently seems that vertex 
corrections applied in the screened Coulomb interaction only 
give the best results. This is corroborated by the fact that this 
quasiparticle dispersion has a better band width and that the 
GqWqT lda introduces little change to the band width. These 
results are in agreement with those of Mahan and Sernelius 52 



6 
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r~zz ~rTT7r~ 

G w r LUA 


r~zz 1 iu i 

GoW L1JA 


Experiment 


(Al) 2.07 


11.5445 


11.6444 


11.1814 


10.60±0.10" 


{Jul) J.ZtS 


4.4044 


A A Q^Q 


A 1 1 OQ 
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(Na) 3.96 


2.9837 


2.9889 


2.7777 


2.65±0.05 c 


(K) 4.96 


1.8625 


1.8579 


1.7044 


1.60±0.05 d 


(Rb) 5.23 


1.6669 


1.6610 


1.5191 


1.70±0.20 e 


(Cs) 5.63 


1.4287 
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FIG. 3: The pair-correlation function — evaluated at r s = 1.0. The 
RPA goes negative for small r, which is a well-known unphysical be- 
haviour of this approximation. Including the vertex makes the Wq DA 
fit the (essentially exact) QMC 51 curve better, but then the p.-c. func- 
tion goes too positive instead. The Wq DA curve does not go to infinity 
when r — > 0, it just goes to a very high value (~ 2000). The horizontal 
lines are there to guide the eye. 
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FIG. 4: Correction to the free-electron quasiparticle dispersion for 
r s = 4.0 — Note the large absolute shift for GoW r LDA . The two 
intersecting straight lines are included to emphasize that none of the 
dispersions fulfill the condition Re \L^ F (£f)] = V xc - 



obtained for a model Hubbard vertex. 



VI. THE CHEMICAL POTENTIAL OF JELLIUM 

To get another indication of whether the large absolute pos- 
itive shift of the quasiparticle dispersion is physical, we com- 
pare with experimental work functions of Al (1 00), (1 1 0) 
and (1 1 1) surfaces (see TablelVTi. We assume that the electron 
density of the surface region, and therefore the electrostatic 
surface-dipole energy barrier, are well described by LDA cal- 
culations including the crystal lattice. The work function, <p, 
will, however, be sensitive to the quasiparticle bulk Fermi 
level, which we use here as a discriminator between self- 
energy approximations in the bulk. Treating the bulk metal 
as jellium, we obtain a shift in the workfunction due to the 



TABLE IV: Occupied band widths of jellium for different r s (eV) - 
evaluated at the self-consistent eigenenergy. 




G W RPA 

G w r LDA 

G W LDA 



FIG. 5: Correction to the free-electron occupied band width — Com- 
parison with experiments. G()Wq' Da has the most significant narrow- 
ing of the band in the relevant density-region. G()W()F LDA is not very 
different from GqW^ pa . Note that a jellium calculation does not in- 
clude the contribution from the crystal lattice potential. 



new chemical potential, 

= <j> LDA + A M = <j> LDA - (Re [E kF (e F )] - V xc ) . 



(16) 



where A/i is the correction due to the shift of the bulk Fermi 
energy for a GW jellium calculation. The LDA workfunc- 
tion is defined as the shift between the vacuum potential, vac , 
and the chemical potential from the LDA surface calculation, 
jU LDA . Since the exact self-energy for jellium must fulfill the 
condition 



fle[£ kF (£ F )]=V X( 



(17) 



we see that the LDA (taken from highly accurate QMC cal- 
culations) corresponds to the exact result if one assumes that 
the bulk is accurately modelled by jellium. Comparing with 
Table [Viand Fig.|6]we see that GoW(f PA is closest to the exact 
result and GoW LDA is slightly further away, while GqWoT lda 
is even worse. 

This leads us to conclude that GoWqT lda is unphysical in 
the sense that a vertex correction derived from a self-energy 
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* Al 


Exp. (eV) c < 


LDA" 


GoW^ 


G () w r LUA 


G < UA 


(100) 


4.41 


-0.14 


0.28 


-1.19 


0.45 


(110) 


4.06 


-0.18 


0.24 


-1.23 


0.41 


(111) 


4.24 


-0.06 


0.36 


-1.11 


0.53 



See reference 
See reference 



TABLE V: The workfunction of aluminium (r s = 2.07) — Compared 
to experiment. The last four columns show the deviation from the 
experimental value. The LDA surface calculation corresponds to 
the exact result if jellium is used to model the bulk fermi energy. 
GoW ( f PA and G % LDA are closer, while GoW r LDA is much worse 



3.5 



3.0 




LDA 



All 1 10) 



Al (100) Al (III) 
•G II WX L 



4.0 4.1 4.2 4.3 4.4 
Experimental Value (eV) 



4.5 



FIG. 6: (Color online) The workfunction of aluminium (r s = 2.07) 
— Compared to experiment. The LDA surface calculation corre- 
sponds to the exact result if jellium is used to model the bulk Fermi 
energy. G W RPA and G W LDA are closer, while G W Q r hDA is much 
worse. The colored lines are there to guide the eye. The black line 
corresponds to perfect agreement with experiment (error bars indi- 
cate experimental uncertainty) 



approximation with a completely local dependence on the 
density (like the LDA) will have pathological features. This is 
most probably due to improper behaviour of the spectral func- 
tion of the screened interaction, as is demonstrated in the final 
section of this paper. 



VII. DISCUSSION AND CONCLUSIONS 

We have presented GqWq DA and GqWqT lda calculations 
for isolated atoms and jellium. We see that GoWqT l wors- 
ens results in all cases compared to the common GqWq FA ap- 
proximation. 

A proper ab initio vertex correction for calculations on an 
arbitrary system should be derived from the starting approxi- 
mation for the self-energy. In this work we have shown that 
in practice, vertex corrections derived from a local density ap- 
proximation to the self-energy (like the LDA) are pathological 
when applied to both the self-energy and the screened inter- 



action. The work function of aluminium was used to confirm 
that the value of the chemical potential in GoWqT lda is far 
from the correct result. 

An indication of why a local correction in both W and E per- 
forms so poorly has been discussed previously by Hindgren 
and Almbladh 58 and investigated in excitonic effects on wide- 
bandgap semiconductors by Marini et al . 59 i 60 . Both types 
of vertex corrections lead to a modified screened interaction 
W = £ v. The spectral function of this, which is required to 
be positive semidefinite for co < and negative semidefinite 
otherwise, is given by 



*q(fl>) 



sgn(o) 



Im [W q (0))] 



sgn(o) Im [£,,(» ]_ 

% |£q(G>)l 



(18) 



so it inherits whatever properties of definiteness the imaginary 
part of the dielectric function has. Now for GoWqT lda this is 
given by, 



Im[e q (a))] = -(v q +K xc )lm [ Zq »] 



(19) 



Since the RPA response function, Xo, will have the correct an- 
alytical properties by construction, this expression will obvi- 
ously change sign whenever K xc - which is strictly negative 
for all densities and a negative constant for jellium - is larger 
in magnitude than v q , which decays as l/q 2 . This will thus 
lead to a spectral function with the wrong sign, which is en- 
tirely unphysical. For jellium, isolated atoms, or any sparse 
enough condensed state, this is guaranteed to happen, because 
K xc — > — °° for low densities. Inspection of the dielectric func- 
tion in G W LDA , 



Im[£ q (<o)] = -- 



t> q Im[ Zq »] 

,x q »| : 



1-^ 



(20) 



illustrates that it cannot suffer from the same pathology. Since 
Eq. d20l > ensures that the static structure factor has the correct 
behavior for both GqWq DA and GqWqT lda , no conclusions 
can be drawn on the reason for the overly positive value of the 
pair-correlation function of jellium when r — > 0, except that it 
must depend on the high k behavior of the denominator. We 
note that none of the calculations have been carried out self- 
consistently; it is possible that the vertices presented here go 
some way to improve self-consistent GW results^. 

One possibility of the failure of the LDA starting point with 
the inclusion of the theoretically consistent vertex is the self- 
interaction error the LDA orbitals carry with them. Any start- 
ing point with an inherent self-interaction error should lead to 
correcting terms in the diagrammatic expansion. It is possible 
that the first-order correction, like GqWoT lda , is not enough 
and higher order corrections must be applied. A vertex derived 
from a second iteration of Hedin's equations does indeed lead 
to further and more complicated diagrams than the equivalent 
vertex from a Hartree starting point. Unfortunately these dia- 
grams are of prohibitive complexity for practical calculations. 

It is still not understood why a correction in W only (in a 
TDDFT-like manner) seems to work as well as it does. There 
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is a similarity here with the way that the Bethe-Salpeter equa- 
tion (BSE) is usually applied for the calculation of optical 
spectra. There too it is well known that, in theory, inclusion 
of a screened interaction in electron-hole excitations should be 
accompanied by an inclusion of the double-exchange term in 
E but this has been proven to worsen results. Recently, Tiago 
and Chelikowsky 62 have used a GqWq DA vertex in conjunc- 
tion with an efficient numerical implementation of the BSE 
for isolated molecules and have shown that the inclusion of a 
TDLDA vertex gives very good results over a wide range of 
structural configurations in excited states. 

For atomic helium and beryllium GqWq DA is very similar 
to the GoW RPA result but slightly worse in neon. While in 
jellium the band width is improved. Hence GqWq DA may be a 
local and easily implementable way to improve quasiparticle 
results in extended systems. 



Overall, vertices based on the local density clearly have 
their limitations, arriving in part from the wavevector indepen- 
dent character of K KC . It should be fruitful to explore vertices 
that incorporate non-local density-dependence and reflect the 
non-local character of the original self-energy operator. 
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